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Abstract 



Using microwave detected, microwave-optical double resonance, we have 
measured the homogeneous linewidths of individual rovibrational transitions 
in the A state of NH 3 , NH 2 D, NHD 2 , and ND 3 . We have used this excited 
state spectroscopic data to characterize the height of the dissociation barrier 
and the mechanisms by which the molecule uses its excess vibrational and 
rotational energies to help overcome this barrier. To interpret the observed 
vibronic widths, a one dimensional, local mode potential has been developed 
along a N-H(D) bond. These calculations suggest the barrier height is roughly 
2100 cm -1 , approximately 1000 cm -1 below the ab initio prediction. The 
observed vibronic dependence of levels containing two or more quanta in V2 
is explained by a Fermi resonance between 2z^ and the N-H(D) stretch. This 
interaction also explains the observed trends due to isotopic substitution. 
The rotational enhancement of the predissociation rates in the NH3 2 1 level 
is dominated by Coriolis coupling while for the same level in ND3, centrifugal 
effects dominate. 
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I. INTRODUCTION 



Although the electronic excited states of most molecules are known to predissociate, the 
exact mechanisms of these processes remain largely undetermined. The difficulty in modeling 
molecular predissociation is that the predissociation is usually caused by a coupling of a 
bound excited electronic state to another electronic state in which the nuclear motion is not 
bound [@ . The rate of dissociation depends upon the properties of both electronic states, 
as well as the nonadiabatic coupling between them. In contrast, the rate of predissociation 
of the ammonia A state is determined by the rate of crossing over, or tunneling through, 
a small barrier on a single potential energy surface ||. While a conical intersection exists 
well beyond this barrier, this intersection effects only the product state distribution and not 
the predissociation rate. Since knowledge of only a single electronic surface is needed to 
understand the A state predissociation, it represents one of the most theoretically tractable 
chemical reactions. 

The nature of the barrier to A state dissociation is well understood. In an early ab initio 
SCF and CI study, the A state was shown to result from excitation to the 4 : a 1 molecular 
orbital (MO) which continuously changed from a predominantly 3s Rydberg orbital at the 
equilibrium conformation to a hydrogen Is atomic orbital at large H2N-H separation ||. A 
small barrier exists along this coordinate due to the change from Rydberg to anti-valence 
character of this orbital as an N-H bond was stretched. The magnitude of this barrier 
was also shown to be a minimum for planar geometry and to increase with out-of-plane 
deformation. The minimum barrier height remains a controversial quantity. In the initial 
ab initio study by Runau et al. |§, the barrier was predicted to lie 4600 cm -1 above the A 
state minimum while a more involved theoretical calculation that allowed the H-N-H angle 
to relax has predicted this barrier to lie at 3200 cm -1 0]. Semi-classical trajectories predict 
it lower by an additional 1200 cm' 1 [||. An accurate determination of this barrier height is 
essential to correctly understand the predissociation mechanisms. Hence, its determination 
was a primary goal of the present work. 

Besides the exhaustive theoretical effort ||-10|, an extensive experimental effort has 



also been directed toward characterizing the barrier height and predissociation mechanisms 
|5|JTT|-p3| . Due to the lack of rotational data as discussed in the previous paper, much 



of the current understanding of the predissociation mechanism(s) has been extracted by 
interpreting the vibronic dependence of the predissociation lifetimes. In particular, Vaida 
and coworkers [20] observed a single fourteen member progression in the absorption spectrum 
of NH 3 and ND 3 in a seeded supersonic jet. Inhomogenous broadening was sufficiently 
reduced such that for each vibronic member of this progression, an average rovibrational 
predissociation rate was extracted. In agreement with the ab initio predictions [|3|, for both 
isotopomers, the predissociation rates were observed to be slowest for the 2 1 level due to 
the strong dependence of the barrier height on out-of-plane deformation. Above this level 
for both species, the predissociation rate was observed be enhanced greatly by increasing z/ 2 
excitation. They postulated that for these higher lying levels anharmonic interactions were 
coupling energy from the bound bending coordinate to the unbound dissociation coordinate, 
thus, circumventing the barrier. Due to the strong anharmonicity of the A state potential, 
only a limited amount of information about the barrier is likely to be revealed through the 
study of these higher lying levels. 
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Deuterium substitution has been shown to reduce the predissociation rate in A state 
ammonia substantially and particularly in the lowest few vibronic levels []12|,|16,20|. Thus, 



determining predissociation lifetimes as a function of the number of hydrogens substituted 
should be a sensitive measure of the barrier height. This effect is attributed to the slower 
tunneling rates of the deuterium compared to hydrogen as well as slower tunneling rates 



due to the lower zero point energy of the N-D stretching modes [24]. If one neglects the 



interaction among the three equivalent N-H(D) bands, the dissociation rate would be the 
sum of tunneling rates through the three equivalent barriers and each subsequent deuterium 
substitution would result in an approximate constant decrease in predissociation rate. Such 
an effect has been observed in the A state action spectra from photodissociation through the 



2 1 level [ 24| . At even greater excitation, coupling of vibrational energy into the dissociation 
coordinate would be expected to enhance the dissociation rate, since it provides a mechanism 
to 'go over the top' of the barrier. Such resonances will depend sensitively on the frequencies 
of the individual modes involved such that a simple isotopic dependence would not be 
expected. The effect of any vibration-rotation interactions on the dissociation rate may also 
lead to a more complex isotopic dependence. 

Since the lowest two vibronic levels could only predissociate by tunneling through this 
barrier, it was anticipated that a better picture of the barrier would be revealed through an 
understanding of the predissociation in the lowest few vibronic levels. Hence, these levels 



have been probed by a host of experimental techniques including absorption [16|, two pho- 



ton fluorescence excitation (2PFE) fl2" |, and rotational resonance Raman excitation profiles 
(RRREP) ||25|| . This work has established a weak predissociation dependence on rotational 
level. This rotational dependence has been quantitatively fit to a model incorporating a 
centrifugal modification to the barrier height |12| ; p5| . However, the much stronger rotational 
dependence observed in the RRREP 's could have also been described by Coriolis forces P5 



These two different mechanisms could not be distinguished since the results of these experi- 
ments represent an average increase among all K sublevels such that distinguishing between 
the two different mechanisms is difficult. 

Fully rotationally resolved spectra have been recorded by optical-optical double reso- 
nance methods P,p2"|. However, these experiments were not able to establish a significant 
predissociation dependence on rotational level. It should be noted that one step of this 
double resonance scheme involves a multiphoton excitation that might be expected to dis- 
tort the observed lineshapes. Thus, due to the experimental uncertainty of the determined 
linewidths in these experiments, a weak predissociation dependence on rotational level might 
have been masked. 

Microwave detected, microwave optical double resonance has recently been used to probe 
the dynamics in the NH3 A state 2° band ||17|| . The rotational dependence observed in this 
level was again quantitatively fit using the centrifugal model In the previous paper, we 



reported the A <— X transition frequencies measured by MODR for the higher lying vibra- 
tional levels of the NH 3 A state as well as for the isotopomers |27|]. This technique has also 
enabled us to measure the homogeneous linewidths for individually resolved Jk rovibronic 
transitions for the entire isotopomeric series. Our observed linewidths are considerably nar- 
rower than previous reports p|,|2"2"||. Thus, it was hoped that these results would enable us 
to discern for the various isotopomers the vibrational and rotational contributions to the A 
state predissociation mechanisms. 
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Assuming that stretching motions are more correctly described by a local mode coordi- 
nate system rather than a normal mode one, a crude one dimensional potential was designed 
along an N-H(D) internal bond coordinate. This potential has been used to demonstrate 
that a 2-1 Fermi resonance between v 2 and V\ quantitatively predicts the observed linewidths 
for levels with more than one quantum in Vi- For the lowest two levels, predissociation could 
only occur by tunneling through this barrier. To reproduce the observed linewidths of the 
vibrationless levels in NH3 and ND3, the barrier height was reduced roughly 1000 cm -1 



relative to the ab initio value. This result agrees with the conclusion of Dixon |28[ who 
also has demonstrated that the ab initio calculations || overestimate the barrier. The pre- 
dissociation efficiency was also shown to be weakly dependent upon rotation. Surprisingly, 
this rotational dependence was found to be isotope dependent. For the hydrogen contain- 
ing isotopomers, the rotational dependence was more adequately described by Coriolis forces 
coupling energy from u 2 and into z/ 3 and z/4, while for ND 3 , this rotational dependence obeyed 
the centrifugal model. This difference arises from the fact that in NH 3 , one quantum of z/ 2 
added to the N-H zero point energy is isoenergetic with scattering states above the barrier 
to dissociation while for ND3, this energy still lies below the top of the barrier. Thus, for 
ND3, the Coriolis contribution to the predissociation rate is reduced such that tunneling 
remains the dominant predissociation pathway. 



II. EXPERIMENTAL 

The recording and analysis of the MODR and RFODR data was described in the previous 
paper and will not be repeated here. The linewidths obtained from such an analysis are given 
in tables |, H, |T|, and [TV]. As suggested by Ziegler |2l| , the observed widths were fit to, 

r = T xexp{T B [J(J+l)-K 2 }+T c K 2 } (1) 

where T is the observed width, r o is a determined constant for a given vibronic band, 
r# reflects the enhancement by rotation in the plane of the molecule, and Tc that of the 
enhancement due to rotation about the symmetry axis. For the asymmetric top species, K 2 
is replaced by (J 2 )/K 2 , calculated from the fitted rotational constants. In these fits, the 
observed widths were weighted by the reciprocal of the squared uncertainty predicted by 
the lineshape fit. Results are given in table |V|. Because the microwave probe transitions are 
Q branch lines of the inversion band, our data is primarily for high K lines. For the mixed 
isotopomer data, we were not able to determine the rotational dependence of the widths 
due to the limited number of transitions observed in these molecules. For the asymmetric 
tops, the ground state Q branch is widely spread, reducing the number of strong transitions 
in the region covered by our microwave equipment. 



III. DISCUSSION 

A. Empirical Local Mode Potential Surface 

The A state absorption spectrum has been modeled as a two symmetric mode problem 
|7|||. One quantum of stretching excitation contains more energy than needed to overcome 



4 



the barrier along a single N-H bond, yet the v\ normal mode coordinate is still strongly 
bound. The dissociation reaction coordinate is a mixture of V\ and the two components 
of the antisymmetric N-H stretching mode, u 3 . This potential topology will lead to strong 
nonlinear interaction between the V\ and u 3 normal modes such that an initial state corre- 
sponding to v\ excitation will rapidly trifurcate into three equivalent wave packets corre- 
sponding to bond-breaking along each of the three equivalent N-H(D) bonds. The depth of 
the exit channels and the sharpness of the motion along the symmetric stretch should lead 



to little amplitude being reflected back. Thus, one should expect only weak recurrence |29 

In work on stable symmetric hydrides, such as the ground state of ammonia, it has 
been established that one can treat the X-H stretching motion equivalently with either a 
local mode or normal mode basis. When the single bond anharmonicity dominates over 
the bond-bond harmonic coupling, the eigenstates are expected to be close to symmetrized 
local mode states. Contrarily, when the harmonic coupling dominates over anharmonicity, 
the eigenstates are closer to normal mode basis functions. For perfect harmonic motion, the 
local mode and normal mode descriptions of the fundamental levels are identical. For the 
ammonia A state, however, the potential along a single N-H bond is highly anharmonic. In 
fact, as will be shown below, the zero point level is the only resonance below the barrier. 
As a result, excitation of one bond mode leads to dissociation in a single vibrational period 
(10 fsec). The time required for exchange of vibrational energy between the N-H bonds is 
longer by approximately f rr j f rr >, which based upon the ab initio harmonic force field |7j is 
predicted to be greater than 50. Thus, the A state of ammonia is the 'ultimate' local mode 
molecule, with bond-bond coupling negligible compared to the rate of dissociation if the 
N-H stretching modes are excited. As a result, we have modeled the dissociation in terms 
of three uncoupled, single bond modes. 

Figure [I] shows a model potential we have used to fit the observed dissociation rate. The 
curve contains information from the ab initio calculations 0] as well as from the experimen- 
tally measured information about the position of the barrier [ 27 1 . These values are repeated 
below. 



r e 


1.06 A 


r b 


1.32 A 


AH 


3200 cm" 1 


f rr (r < r b ) 


2.56 x lO^m" 1 /! 2 


f„{r b <r <R) 


-3.36 x Wcm^/A 2 



The curve consists of a pair of parabolas whose extrema and curvature are based upon the ab 
initio potential Q, but whose barrier height can be adjusted. Beyond some large separation, 
R S> r e , the potential was taken to be -10000 cm -1 correlating with the energy difference 
of the NH 3 A state and the ground states of the photofragments, NH 2 and H |§. The real 
eigenstates, ?/>(E) were calculated at each energy by numerical integration and normalized 
such that fip(E)ip(E')dx = 2n5(E - E'). 

We begin our discussion with the barrier height picked to agree with the most recent ab 
initio prediction of 3200 cm -1 [|]]. In order to mimic what we would expect for V\ absorption 
activity, we calculated the Franck-Condon overlap of ^(E) with a Gaussian picked to model 
the zero point level of the ground electronic state. The resulting squared overlap versus 
energy is presented in figure |2|. We observe a sharp resonance at 1450 cm -1 with a HWHM of 
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0.16 cm -1 . This resonance corresponds to the zero-point level below the dissociation barrier. 
If this represented the vibrationless level of the A state, one predicts a homogeneous width 
of three times this size since there are three equivalent hydrogen atoms that could tunnel 
out; thus, the lifetime broadening will be three times the lifetime broadening of each N-H 
bond. This prediction is much smaller than the 17 cm' 1 HWHM measured for the NH3 
origin. 

While the features of our crude potential have a clear physical interpretation, this po- 
tential overestimates the thickness of the barrier relative to the ab initio potential J| . As a 
consistency check, using a spline interpolation, we numerically solved the one-dimensional 
Schroedinger equation for this potential. Franck-Condon overlap with the ground state 
wavefunction resulted in a sharp resonance at 1460 cm" 1 with a width of 0.73 cm -1 in 
qualitative agreement with the results from our more crude approach. Thus, we justify our 
simple method and conclude that the ab initio calculation has overestimated the barrier 
height. Dixon reached the same conclusion based upon his two-dimensional wavepacket 
calculations ||. 

In order to reproduce the observed 17.1 cm -1 width of the NH3 zero point transition, the 
barrier height in our model potential was reduced to 2075 cm" 1 . Plotted in figure ^| is the 
Franck-Condon overlap as a function of energy above T e for this corrected potential. The 
zero point level lies 1350 cm -1 above T e and has a HWHM of 5.7 cm" 1 per hydrogen in good 
agreement with the experimental measurement. Thus, in order to model the experimentally 
measured width of the NH3 zero point level, we need to lower the ab initio barrier height 0] 
by more than 1000 cm" 1 . 

As shown in figure^, an extremely broad quasi-resonance above the barrier is also present. 
This quasi-resonance represents what may be called one quantum in the N-H local mode. 
Eigenstates in this energy region have their inner turning points at bond lengths slightly 
less than r e , and thus have significant overlap with the ground state Gaussian. Our crude 
calculation predicts this resonance lies 3370 cm -1 above T e (2020 cm -1 above the zero point 
level) and has a 2850 cm -1 linewidth. Because of its broad width, the peak absorption is 
expected to be small relative to the zero point level, making this level very difficult to detect. 
Ashfold et al. || have measured a 500 cm" 1 wide resonance centered at 2300 cm" 1 above 
the zero point level that they have assigned to l 1 . As discussed in our previous paper, we 
have reassigned this band to 4 2 . This level is very likely to be strongly coupled to the N-H 
stretching modes through Fermi coupling leading to the observed rapid dissociation rate. 

Using the same potential, we also calculated the Franck-Condon overlap for an N-D 
bond. Results from the calculation are presented in figure f|. The zero point level is found at 
1035 cm" 1 with a width of 0.16 cm -1 per deuterium. The ND 3 zero point level is predicted 
to have a homogeneous linewidth of 0.48 cm" 1 comparing favorably with the experimental 
measurement of 0.4 cm -1 . We calculate the N-D second resonance to have a homogeneous 
linewidth of 1360 cm -1 and to lie about 2750 cm -1 above T e or about 1720 cm -1 above 
the zero point level. A broad feature was observed at 1790 cm -1 in the ND 3 dispersed 
fluorescence spectrum that has been assigned to the l 1 band. As for NH 3 , we reassign this 
feature to the 4 2 band. 

A quantitative calculation to test the accuracy of this local mode approximation is clearly 
desirable. Unfortunately, this calculation would require at a minimum ab initio calculations 
along the three stretching coordinates, and then a 3-D wavepacket propagation of the ground 
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state Gaussian on this potential. If such a calculation dramatically reduces the predicted 
resonance widths of the l 1 level, then we would withdraw our reassignment of the peaks 
Ashfold et al. |5j assign to l 1 . However, since the basic vibronic and isotopic trends are 
predicted by this crude potential, we believe our empirical dissociation potential contains 
the essential features of the true potential. 



B. Predissociation dependence upon vibrational level 

The increasing dissociation rate with v 2 excitation above 2 1 suggests that anharmonic 
interactions are coupling two quanta of vibrational energy out of the bending normal mode 



into one of the three N-H(D) bond coordinates lying along the dissociation coordinate |20 
The lowest order vibrational coupling that could be responsible is a Fermi resonance between 
v 2 and V\ caused by the khhs term in an harmonic approximation. We estimate the value of 
this coupling constant by assuming the molecule vibrates out-of-plane with constant bond 
length. This approximation is justified by the ab initio potential Q. 

Classically, the above coupling leads to a forcing of the stretching mode at two times the 
frequency of the bending mode. If there is a 2:1 resonance of the harmonic frequencies, one 
would expect an efficient transfer of energy inside some resonance zone. The experimentally 
determined force field given in the previous paper does not predict an exact resonance p7| , 



but a small amount of average energy transfer would be expected due to non-resonant cou- 
pling. Further, given the rapid dissociation of the N-H(D) stretching modes, their frequencies 
are poorly defined, making the resonance condition much less restrictive. A 2:1 Fermi reso- 
nance implies the transfer of two quanta in v 2 into one quantum in V\\ anharmonic transfer 
of one quantum of v 2 is rigorously symmetry forbidden. We see that this mechanism oper- 
ates only for levels containing two or more quanta in v 2 and the rapid increase in linewidth 
observed above the 2 1 level is predicted. 

We use our empirical potential to estimate the expected magnitude of the Fermi reso- 
nance. We use as a zero order model, the coupling of two harmonic oscillators and by second 
order perturbation theory, we find that the average vibrational energy in the stretching co- 
ordinate increases by an amount given by, 



(AE S ) = u s 



( khhs \ 

UV2J 



u 2 (u 2 -l) (2z/ 2 + l) 2 (i/ 2 + l)(z/ 2 + 2) 



+ 



(co s - 2uo 2 ) 2 to 2 (uj s + 2uj 2 ) 2 



(2) 



where uj s is the stretching frequency, uj 2 is the bending frequency, v 2 is the number of quanta 
in the out-of-plane bend, and k hhs is the coupling constant. Using the results from Ziegler 



2q| , the effect of the Fermi resonance on the transition linewidth is 

(ae s )- 



AT = T x exp 



(3) 



where %uo 2 = 1063 cm -1 is the imaginary frequency at the top of the barrier. Putting in 
values estimated from the ab initio force field [|7| gives an increase in energy that is far too 
small to explain the observed vibrational dependence of the widths. 

In order to account for the highly anharmonic motion along the N-H(D) stretching 
coordinates, we use our above estimate for this potential curve. States with two quanta in 
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the out-of-plane bending mode couple directly to scattering states above the dissociation 
barrier. We use Fermi's Golden rule to estimate the contribution of the coupling to the 
dissociation rate, 

k diss = f (4) 

where (Vo) is the coupling matrix, E is the energy relative to the bottom of the exit channel, 

and the \j2m/E factor arises due to the density of states and box normalization of the 
wavef unctions. Assuming the simplest resonance, Fermi coupling between the bend and the 
stretch is predicted by 

(V c ) = \frhh{v2 - 2\h 2 \u 2 )(E + AE 2 \(r - r )\^ ) 2 (5) 

where f r hh is the cubic bend stretch force constant, h is the out-of-plane bend internal 
coordinate, and (E Q + AE 2 \(r — r )\^ ) is the coupling matrix element. \^ ) is the zero 
point resonance normalized to unity inside the well and (E a + AE 2 \ is a continuum function 
normalized to 2%sin(kr) in the asymptotic region (r ^> r ). We use the harmonic approxi- 
mation for the out-of-plane bend and values for fhh from our experimentally estimated force 
field to determine the expectation value for h 2 . This value was given by 

(y 2 - 2\h 2 \u 2 ) = 7.614 x 10~ 3 y/u 2 (u 2 - 1)A 2 (6) 

Given in figures [5] and ^ is the squared coupling matrix element as a function of the energy 
in the N-H(D) bond. Collecting the above expressions along with the ab initio value for 
f r hh of 2.66 aJ/A 3 allowed estimation of the predissociation rates and thus the expected 
transition linewidths. These calculated values are presented in table [VT| along with our 
measured values. 

The predicted dissociation rates are about a factor of three too large for the 2 2 levels, and 
increase too rapidly above that. Roughly, the isotopic dependence is correctly predicted. 
One explanation for the predicted dissociation being much faster than observed is that the 
ab initio value of f r hh may be too large. Even more significant may be our neglect of 
the barrier dependence on out-of-plane angle. Increasing the barrier height will reduce the 
coupling matrix element among the zero point level and the scattering states and decrease the 
predicted Fermi resonance contribution to the dissociation rate. Consequently, we conclude 
that our overestimate of the Fermi contribution to the dissociation rate is caused by the 
neglect of the barrier dependence on out-of-plane deformation. Proper treatment of this 
problem, including both Fermi resonance and the effect of barrier changes with out-of-plane 
angle, requires two-dimensional wavepacket calculations that are beyond our computational 
resources. 

This resonance calculation does rationalize the larger than expected dissociation rate for 
the NHD2 isotopomer. Even though there is great uncertainty as to the frequency of the 
stretching levels, two quanta of v 2 for NHD2 would be expected to be in closer resonance 
with the N-D stretch than for two quanta in any of the other isotopomers, leading to an 
enhanced dissociation rate. Contrary to the mechanism controlling the dynamics in the 
tunneling region, this mechanism predicts that the N-D bond dissociates faster than the 
N-H bonds for these higher lying levels. 
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C. Predissociation dependence on rotational level 



The predissociation dependence on rotational level can arise from two independent 
sources. In the first mechanism, as extensively discussed by Ashfold et al. [0,0] an d Ziegler 



2q| , rotation about either the b or c axes will promote an N-H(D) bond lengthening result- 
ing in an effective barrier reduction. In the second mechanism, b axis rotation causes the v% 
mode to be coupled to the unbound u 3 and 1/4 modes through Coriolis forces and, hence, in- 
creasing the probability of barrier penetration. If centrifugal effects are the dominant source 
of rotational enhancement, the ratio of b to c axis effects should be constant with vibronic 
excitation. Contrarily, the Coriolis coupling matrix elements would scale with z/ 2 . Hence, 
the two mechanisms can be distinguished by measuring the ratio of b to c axis effects as a 
function of z^. 

We estimate this ratio allowing only centrifugal effects to contribute. Assuming the 
dissociating H(D) atom lies on the y axis, the effective potential is given by 

Y t l >\ T T / >\ \J{J+1)-K 2 } K 2 

v(r) = v ° {r) + 44(?) + mr) (7) 

where V (r) is the rotationless potential and 1/2[J( J+ 1) — K 2 ] and K 2 are the expectation 
values for x and z axis rotation respectively. Using the ab initio values for r e and r b at 
1.055 A and 1.32 A respectively, the rotational correction to the barrier height is estimated 
to be 

AV = 1.25[J(J+ 1) - K 2 } +0.73K 2 . (8) 

Thus, the centrifugal model predicts a B axis effect to C axis effect ratio of 1.75. We also 
note that for the centrifugal model this ratio is independent of the barrier shape and only 
weakly dependent upon barrier position. For instance, as was varied from r e to oo, this 
ratio only changed from 2 to 1. 

From our NH 3 data, the ratio of b axis to c axis effects was found to be 2.2 ± 1. for the 
zero point level and 15. ± 2. for the 2 1 level. We predict that centrifugal effects dominate 
in the vibrationless level while Coriolis forces dominate in the 2 1 level. This change in 
dominant mechanism can be explained utilizing our empirical N-H local mode potential. 
This potential predicts the zero point energy of an N-H bond mode to be 1350 cm -1 while 
the barrier height was determined to be 2075 cm^ 1 . The 2 1 level lying roughly 900 cm -1 
above this level is expected to lie well above the barrier and therefore be Coriolis coupled 
to a manifold of continuum states. 

Using a procedure similar to the one we used to estimate the magnitude of the Fermi 
resonance, we estimate the rotational dependence of the linewidths in the NH3 2 1 level due 
to Coriolis forces. In this case, the coupling matrix elements are given by 



(Vc) = -I-^J{J + l)-K 2 Q h (Q s P b - Q b P s ) (9) 

where I xx is the moment of inertia about the x axis, Q b is the stretch-bend Coriolis coupling 
constant which we set to its ground state value of 0.752, Q s and Qb are the normal mode 
coordinates for the stretch and bend, and P s and Pb are the conjugate momenta to these 
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normal coordinates. By substituting equation |9] into equation |, we can estimate the Coriolis 
enhancement to the predissociation rate and the expected linewidth enhancement. The 
Coriolis contribution to the linewidth was calculated to be 

r o = 0.n[j(j+i)-^ 2 ] (io) 

For the NH3 2 1 level, this predicted value is in excellent agreement with the observed b axis 
rotational dependence of 0.10(4) cm -1 . 

The ND 3 results are more poorly determined than our NH 3 results. We still can compare 
our determined ratio of b axis to c axis enhancement for the two isotopomers. In contrast 
to NH 3 , this ratio for the ND 3 2 1 level was determined to be 2.1(7) more in line with 
the centrifugal model. Surprisingly, for the two isotopomers, the dominant mechanism for 
rotational enhancement appears to be different. 

To understand this result we again turn to our empirical N-H local mode potential. This 
potential predicts the zero point energy in the N-D stretch to be 1035 cm -1 . Addition of 
one quantum in with an energy of approximately 670 cm -1 still lies well below the barrier. 
Since continuum states below the top of the barrier have small amplitude near r e , the Coriolis 
coupling matrix element should be much smaller for ND 3 than for NH 3 . Therefore, changes 
in the zero point energies and mode frequencies due to isotopic substitution are shown to 
have a dramatic effect on the predissociation dependence on rotational level. 



IV. CONCLUSION 

We have measured the homogeneous linewidths of individual rovibrational transitions 
using MODR. Our measured linewidths are considerably narrower than previous determi- 
nations P,|2"2"[ . Since the observed widths were much larger than the laser bandwidth or the 



Doppler width, they represent a direct measure of the predissociation lifetime. Thus, we 
have used these linewidths to more fully interpret the A state predissociation mechanisms. 

Using the available experimental and theoretical data, we have developed a one- 
dimensional, local mode potential along a N-H(D) bond. This potential has been used 
to qualitatively predict the vibronic and isotopic dependence of the predissociation rates. 
In order to simulate the observed widths in the zero point levels, the barrier height had to 
be lowered to 2075 cm -1 , about 1000 cm -1 below the ab initio prediction ||. Dissociation 
rates from the 2° and 2 1 levels could only proceed by tunneling through this barrier while 
dissociation from the higher lying levels was described by a 2:1 Fermi resonance of ^2 di- 
rectly to the manifold of N-H(D) continuum states lying above the barrier. This model also 
predicts that photdissociation of NHD 2 should proceed primarily through the NHD + D 
pathway. 

Dissociation of the two lowest vibronic bands was found to have a weak dependence 
upon rotation level. Rotationally enhanced dissociation from the zero point level in both 
NH 3 and ND 3 is caused by the lowering of the effective barrier height by centrifugal forces . 
The rotational enhancement of the dissociation rate from the NH 3 2 1 level, which is predicted 
to lie above the potential barrier along the N-H stretch, was dominated by perpendicular 
Coriolis coupling to the manifold of continuum states lying above the barrier. The same 
level in ND 3 lies below the barrier to dissociation and the rotational dependence of the 



dissociation rates still fits the centrifugal model • 
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TABLES 

TABLE I. Linewidths of the NH 3 A <- X state MODR transitions. 



2° band. 



Transition 


HWHM (cm x ) 


la 


Q 2 (2) 


16.83 


0.06 


R 2 (2) 


19.38 


0.21 


Qs(3) 


17.90 


0.05 


R 3 (3) 


19.37 


0.22 


Q 4 (4) 


18.81 


0.11 


Q 4 (5) 


18.11 


0.25 


R 4 (5) 


24.20 


3.60 


Qs(5) 


19.19 


0.13 


R 5 (5) 


25.98 


0.99 




1 n a q 
ly.4o 


n nc 
U.Uo 


Re(6) 


33.31 


1.52 


P 6 (7) 


16.27 


1.57 


Qe(7) 


20.79 


0.14 


R 6 (7) 


26.48 


1.37 


Q 7 (7) 


21.04 


0.14 


R 7 (7) 


23.70 


1.88 


Qs(8) 


21.00 


0.11 


R 8 (8) 


26.03 


2.35 


Qg(9) 


22.74 


0.44 


2 1 band. 


Transition 


HWHM (cm- 1 ) 


la 


Qi(i) 


15.03 


0.07 


Ri(l) 


14.82 


0.15 


P 2 (3) 


12.73 


0.16 


Q 2 (3) 


15.39 


0.05 


R 2 (3) 


15.99 


0.08 


Qs(3) 


15.07 


0.02 


R 3 (3) 


15.48 


0.09 


P 2 (5) 


13.88 


0.33 


Q 2 (5) 


15.56 


0.45 


R 2 (5) 


17.24 


0.27 


P 4 (5) 


13.64 


0.14 


Q 4 (5) 


19.65 


0.16 


R 4 (5) 


16.35 


0.06 


Q 4 (5) 


14.75 


0.43 


R 4 (5) 


17.68 


0.31 


P 3 (6) 


18.27 


0.35 


Q 3 (6) 


12.10 


0.81 


R 3 (6) 


16.75 


0.08 
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P 5 (6) 


18.98 


0.54 


Q 5 (6) 


11.09 


0.54 


R 5 (6) 


17.14 


0.08 


Re(7) 


21.28 


0.54 


Qt(8) 


17.93 


0.15 


Qs(8) 


16.04 


0.03 


R 8 (8) 


20.10 


0.79 


P 8 (9) 


46.17 


6.61 


R 8 (9) 


23.22 


2.33 


Q 9 (9) 


14.64 


0.09 


R 9 (9) 


33.74 


5.11 


2 2 band. 


Transition 


HWHM (cm x ) 


la 


Qs(3) 


22.08 


0.27 


R 3 (3) 


26.02 


1.36 


Qe(6) 


20.09 


0.16 


R 6 (6) 


21.92 


4.07 


2 3 band. 


Transition 


HWHM (cm- 1 ) 


la 


Qi(i) 


42.37 


0.58 


Q 2 (2) 


45.30 


0.64 


Qs(3) 


45.33 


0.30 


Qs(4) 


47.32 


0.72 


Q 4 (4) 


53.33 


0.35 


Q 3 (5) 


48.36 


0.23 


Q 5 (5) 


49.66 


0.53 


Qs(6) 


47.94 


0.71 


Qe(6) 


46.91 


0.33 


Qe(7) 


48.53 


0.22 


Qr(7) 


44.54 


0.24 


Qe(8) 


45.41 


0.61 
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TABLE II. Linewidths of the NH 2 D A <- X state MODR transitions. 



2° band. 



Transition 


HWHM (cm- 1 ) 


la 


3o3 3i3 


17.46 


0.10 


^23 <~ 3i3 


17.46 


0.10 


4i4 <- 4 04 


17.41 


0.10 


5 14 <- 4 04 


17.41 


0.10 


5l5 5o5 


18.57 


0.18 


615 5o5 


18.57 


0.18 


606 7i6 


22.85 


0.46 


^26 7i6 


22.85 


0.46 


7o7 <— 817 


18.97 


0.18 


827 817 


18.97 


0.18 


927 817 


19.47 


1.32 


2 1 band. 


Transition 


HWHM (cm- 1 ) 


la 


3l3 3o3 


11.74 


0.04 


4l3 3o3 


12.27 


0.26 


4 4 <- 4u 


12.21 


0.14 


5 2 4 <- 4u 


12.21 


0.14 


5o5 5i5 


12.43 


0.09 


625 <- 5i 5 


20.33 


1.53 


6l6 726 


10.35 


1.11 


7l6 726 


13.54 


0.14 


836 *— 726 


11.83 


0.73 


7l7 827 


11.64 


0.94 


817 <— 827 


13.74 


0.08 


937 <- 827 


15.27 


0.74 


2 2 band. 


Transition 


HWHM (cm" 1 ) 


la 


3o3 3i3 


21.30 


0.12 


423 3i3 


21.30 


0.12 


4i4 <— 4o4 


19.60 


0.28 


5l5 5o5 


21.85 


0.13 


615 <— 5o5 


21.85 


0.13 


726 7i6 


20.48 


0.78 


7o7 <— 817 


7.83 


3.43 


827 817 


21.62 


0.40 


927 817 


42.11 


7.39 
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TABLE III. Linewidths of the NHD 2 A <- X state MODR transitions. 



2° band. 

Transition HWHM (cm" 1 ) To~ 

2 02 <- 2i2 1024 Oil 

3 22 <- 2i2 11.12 0.35 

3i 3 <- 3 03 9.41 0.17 

4i 3 <- 3 03 6.58 0.83 

4 W <- 5 24 8.18 0.99 

5i4 <- 5 24 10.30 0.12 

6 34 <- 5 24 10.09 0.45 

2 1 band. 

Transition HWHM (cm -1 ) lcT~ 

2i2 <- 2 02 045 O04 

312 <- 2 02 5.95 0.14 
3 03 <- 3i 3 6.29 0.05 
4 23 <- 3i 3 6.35 0.28 
4 04 <- 5 M 6.66 0.21 
5 24 <- 5i4 6.54 0.03 
6 24 <- 5i4 7.13 0.10 

2 2 band. 

Transition HWHM (cm -1 ) To~ 

2 02 <- 2i2 17.88 018 

3 22 <- 2i2 17.88 0.18 

31 3 <- 3 03 18.98 0.28 
4i 3 <- 3 03 18.98 0.28 
4 W <- 5 24 18.99 0.09 
5i4 <- 5 24 18.99 0.09 
6 34 <- 5 24 18.99 0.09 
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TABLE IV. Linewidths of the ND 3 A <- X state MODR transitions. 

2° band. 



Transition 


HWHM (cm- 1 ) 


la 


Qs(3) 


2.241 


0.099 


Q 4 (4) 


1.927 


0.051 


Qs(5) 


1.916 


0.096 


Q 5 (5) 


2.024 


0.036 


Qe(6) 


1.946 


0.075 


Qe(7) 


1.364 


0.114 


Qr(7) 


2.085 


0.038 


Qe(8) 


1.790 


0.056 


Qr(8) 


2.193 


0.053 


Qs(8) 


2.075 


0.026 


Qe(9) 


1.581 


0.355 


Q 9 (9) 


1.524 


0.062 


Q 9 (10) 


1.887 


0.070 


2 1 band. 


Transition 


TT"\T5 7"TTT\ IX { ~. — 1\ 

HWHM (cm J 


-t 

1(7 


Qi(i) 


0.782 


0.044 


Q2(2) 


0.781 


0.030 




0.664 


0.016 


R 3 (3) 


0.742 


0.029 


Q 3 (4) 


0.749 


0.025 


Q 4 (4) 


0.765 


0.008 


R 4 (4) 


0.743 


0.059 


Q 3 (5) 


0.832 


0.053 


Q 5 (5) 


0.936 


0.033 


Qs(6) 


1.080 


0.053 


Q 5 (6) 


1.050 


0.036 


Qe(6) 


1.518 


0.019 


Qe(7) 


0.958 


0.017 


Qr(7) 


1.049 


0.014 


Qe(8) 


1.344 


0.021 


Q 7 (8) 


1.106 


0.050 


Qs(8) 


1.039 


0.012 


Q 3 (9) 


1.540 


0.137 


Qe(9) 


1.923 


0.127 


Q 9 (9) 


1.281 


0.031 


Q 9 (10) 


1.611 


0.029 


Qg(ll) 


1.849 


0.058 


Qio(H) 


1.593 


0.075 


Qn(12) 


1.598 


0.058 


Qi 2 (14) 


1.534 


0.086 


2 2 band. 
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Transition 






HWH1V1 (cm 


) 


1 

1(7 


Qs(3) 






7.226 




0.010 


Q3(4) 






1 1 r AO 

11.592 




0.328 


04(4) 






T n C\ 

7.ofc>0 




0.198 








9.423 




0.463 


Qs(5) 






9.941 




0.172 


Qe(6) 






7.640 




0.148 


Qe(7) 






6.814 




0.505 


Qt(8) 






8.789 




0.221 


Q 9 (9) 






8.352 




0.704 


TABLE V. 


Constants obtained from a 


fit of the observed rovibronic linewidths to equation p]. 




^2 


r D 




r B 


r c 


TVTTT 

JNH 3 





17.0(4) 




0.006(4) 


0.0028(9) 




1 


14.6(3) 




0.007(2) 


0.0004(5) 




2 


22(1) 




0.02(1) 


-0.005(2) 




3 


47(4) 




0.001(4) 


-0.001(2) 







17.0(9) 










1 


11.5(3) 










2 


21(1) 








IN rlL)) 





10(1) 










1 


6.4(3) 










2 


17.7(5) 








ND 3 





2.1(1) 




-0.0048(29) 


-0.0016(11) 




1 


0.74(5) 




0.0058(31) 


0.0038(15) 




2 


7.7(8) 




0.065(11) 


-0.0021(47) 
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TABLE VI. Linewidths calculated for 2:1 Fermi coupling between the out-of-plane bend and 
the stretching vibration. 



NHs 



V-2 

2 
3 
4 



T(calc) 



r(obs) 



66.8 
160.4 
244.1 



22 
47 
75 



NH 2 D 

2 79.0 21 

3 191.5 52 

4 313.3 78 
5 451.5 78 

NHD 2 

2 77.7 18 

3 207.5 36 

4 343.2 54 
5 474.6 59_ 

ND 3 

2 37.3 7.6 

3 98.5 

4 162.5 

5 267.2 17.5+ 
tReference Wi 
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FIGURES 



FIG. 1. Model one-dimensional energy surface of the A state of ammonia along the dissociative 
N-H coordinate. Curve is composed of two parabolas with positions and radius of curvature 
corresponding to the calculated values. 

FIG. 2. Predicted Franck-Condon overlap for N-H from ground state wavefunction. Upper 
state potential has a 3200 cm -1 high barrier to dissociation. First excited state lies above the 
barrier and lies 2480 cm -1 above the zero point level and has <1% of the peak height. 

FIG. 3. Predicted Franck-Condon overlap for N-H from ground state wavefunction as a function 
of energy above T e . Upper state potential has a 2075 cm -1 high barrier to dissociation. The first 
excited state lies appoximately 2000 cm -1 above the zero point level and is approximately 950 cm -1 
wide. 

FIG. 4. Predicted Franck-Condon overlap for N-D from ground state wavefunction as a function 
of energy above T e . Upper state potential has a 2075 cm -1 high barrier to dissociation. The furst 
excited state lies appoximately 1700 cm" 1 above the zero point level and is approximately 450 cm -1 
wide. 

FIG. 5. The squared bend-stretch coupling matrix element for the bend and the stretch versus 
the energy in a single N-H stretch including the zero point energy. 

FIG. 6. The squared bend-stretch coupling matrix element for the bend and the stretch versus 
the energy in a single N-D stretch including the zero point energy. 
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Figurel 
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